Analysis of Dynamic Multiplicity Fluctuations at 
PHOBOS 

Zhengwei Chai 2 for the PHOBOS Collaboration 
B B Back 1 , M D Baker 2 , M Ballintijn 4 , D S Barton 2 , R R Betts 6 , A 
A Bickley 7 , R Bindel 7 , A Budzanowski 3 , W Busza 4 , A Carroll 2 , Z 
Chai 2 , M P Decowski 4 , E Garcia 6 , N George 1 ' 2 , K Gulbrandsen 4 , S 
Gushue 2 , C Halliwell 6 , J Hamblen 8 , G A Heintzelman 2 , C 
Henderson 4 , D J Hofman 6 , R S Hollis 6 , R Holyhski 3 , B Holzman 2 , A 
Iordanova 6 , E Johnson 8 , J L Kane 4 , J Katzy 4,6 , N Khan 8 , W 
Kucewicz 6 , P Kulinich 4 , C M Kuo 5 , W T Lin 5 , S Manly 8 , D 
McLeod 6 , A C Mignerey 7 , R Nouicer 6 , A Olszewski 3 , R Pak 2 , I C 
Park 8 , H Pernegger 4 , C Reed 4 , L P Remsberg 2 , M Reuter 6 , C 
Roland 4 , G Roland 4 , L Rosenberg 4 , J Sagerer 6 , P Sarin 4 , P Sawicki 3 , 
W Skulski 8 , P Steinberg 2 , G S F Stephans 4 , A Sukhanov 2 , J -L 
Tang 5 , A Trzupek 3 , C Vale 4 , G J van Nieuwenhuizen 4 , R Verdier 4 , F 
L H Wolfs 8 , B Wosiek 3 , K Wozniak 3 , A H Wuosmaa 1 , B Wyslouch 4 

1 Argonne National Laboratory, Argonne, IL 60439-4843, USA 

2 Brookhaven National Laboratory, Upton, NY 11973-5000, USA 

3 Institute of Nuclear Physics, Krakow, Poland, 

4 Massachusetts Institute of Technology, Cambridge, MA 02139-4307, USA 

5 National Central University, Chung-Li, Taiwan 

6 University of Illinois at Chicago, Chicago, IL 60607-7059, USA 

7 University of Maryland, College Park, MD 20742, USA 

8 University of Rochester, Rochester, NY 14627, USA 



Abstract. This paper presents the analysis of the dynamic fluctuations in the inclusive 
charged particle multiplicity measured by PHOBOS for Au+Au collisions at ^/sjvjv = 200GeV 
within the pseudo-rapidity range of —3 < r\ < 3. First the definition of the fluctuations 
observables used in this analysis is presented, together with the discussion of their physics 
meaning. Then the procedure for the extraction of dynamic fluctuations is described. 
Some preliminary results are included to illustrate the correlation features of the fluctuation 
observable. New dynamic fluctuations results will be available in a later publication. 



One of the important goals of the experiments at the Relativistic Heavy Ion Collider (RHIC) 
is to search for the evidence of the quark-gluon plasma(QGP) /hadron gas(HG) phase transition. 
Event-by-event fluctuations might be a useful signature of phase transition evidence, since they 
could be significantly altered if a phase transition occurs shortly after the collision P0| due to 
the the large difference between the degrees of freedom of the two phases. Since fluctuation 
observables are intrinsically related to particle correlations, the study of fluctuations could also 
provide helpful insights on the mechanisms of particle production in heavy ion collisions. 

In this analysis, we measure the inclusive charged particle multiplicity fluctuations as a 



function of the separation and bin size in pseudo-rapidity^) space in Au+Au collisions with an 
energy of ^snn = 200GeV. Since the event-by-event particle multiplicity is mainly determined 
by the the number of participating nucleons (N par t) in the collision, a significant amount of 
the multiplicity fluctuations is due to the N par t variation, which is just a geometry effect. To 
study the dynamic fluctuations which are related to intrinsic particle production mechanisms, 
the Np ar t fluctuations must be removed from the measurement of multiplicity fluctuations. 

For the study of multiplicity fluctuations, we use an event-by-event observable C, which is 
defined as 

c= N ^~ N \ (1) 

where N\ and N 2 are the multiplicities in a pair of r\ bins with the same bin size and symmetric 
with respect to 77 = 0. The width of the C distributions (<r(C)) is used as our fluctuations 
observable. Because the difference of the two multiplicities is used in the definition of C and 
the two multiplicities change in the same direction when N par t varies from event to event, 
the Np ar t fluctuations are suppressed in the measured cr(C). Also because of event-by-event 
normalization factor y/N[ + N 2 in the denominator, <r(C) = 1 for independent particles. Thus 
non-1 cr(C) indicates non-0 dynamic fluctuations in addition to the statistical fluctuations. 
We can decompose the multiplicity fluctuations c 2 (C) into two parts: statistical fluctuations 
( a stat) an d dynamic fluctuations (cr^ ). While the statistical fluctuations are due to the finite 
multiplicity, the dynamic fluctuations are related to the intrinsic correlations in the particle 
production. The long range (particles from different bins) correlations and the short range 
(particles usually in the same bin) have different effect in cr 2 (C) and different signs of the 
dynamic fluctuations. 

To understand what effect the long range correlations have on cr 2 (C), we can use the 
approximate expression for a 2 {C): c 2 (C) a ( N ^)+ a ^ 1 ^)-^P^( N ^) a ( N ^) ~ 1 — ^ where p is 
the correlation coefficient between N\ and N2, assuming the multiplicity N follows Poisson 
distribution (cr(N) = \/N). With such approximation, <j\ yn ^ —p, positive long range 
correlations lead to negative dynamic fluctuations and reduce cr 2 (C) from 1. The short range 
correlations can be treated as a cluster effect and related to the UA5 cluster multiplicity study [2J. 
Assuming the particles are created in cluster with multiplicity k, and the multiplicity in the 
investigated bin is increased by a factor of k, C = is increased to y/kC = 7==== and 

<t(C) is changed to Vka(C) correspondingly. We then have <7 2 (C) = <?1tat + a dyn = & which is 
a measure of cluster multiplicity, if all the associated particles in the cluster are included in one 
bin. In practice not all the associated particles in cluster can be included in the same bin, and 
<t 2 (C) can be regarded as an "effective" cluster multiplicity. Thus the short range correlations 
lead to positive dynamic fluctuations and increase <J 2 (C) from 1. There is also the case where the 
short range correlations show long range effect in a 2 (C), if the correlated particles are selected 
into different but close bins. 

For cr(C) calculated from reconstructed event multiplicity, it also includes additional detector 
effects. Since both the C distribution of the particles from the model event generator and the C 
distribution of the reconstructed multiplicity after detector response simulation are very close 
to Gaussian distributions, we may treat the detector effects as Gaussian response smearing and 
write the reconstructed multiplicity fluctuations <r 2 (C) as the sum of three parts: 

a 2 (C) = a 2 tat + a\ yn + aj et , (2) 

where the a 2 det represents the detector effects. Using the detector simulation to estimate o\ eV 
the dynamic fluctuations can be extracted from the following simple equation: 

Vdyn = ° 2 (C) - Vstat ~ a def ( 3 ) 
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Figure 1. <r(C) versus r\ with fixed r\ bin size 0.5. All values of <r{C) were calculated using 
50% azimuthal acceptance. The estimated systematic error of 5% is not shown in the plot. 



This analysis used data taken by the PHOBOS detectors during the runs providing Au+Au 
collisions at ^snn = 200GeV. Only those data taken with a zero magnetic field setting and 
minimum bias triggers were used. The PHOBOS octagonal multiplicity detector was used to 
measure the charged particle multiplicity. It consists of 92 silicon sensors which surround the 
beam pipe and cover the pseudo-rapidity range of \rj\ < 3. The active elements of the octagonal 
detector are constructed of highly segmented silicon pads which yield signals at the passage of 
charged particles. Within this detector, every second azimuthal octant has a window for another 
detector (spectrometer or vertex detector), which leads to incomplete azimuthal coverage. To 
reduce the systematic uncertainties in the fluctuations results, these nonuniform octants were 
excluded in the analysis of the reconstructed multiplicities from both the data and simulations, 
which left a total nominal azimuthal coverage of 50%. 

In each rj bin, the sum of the angle-normalized energy deposits (dE/dx) was used to 
estimate the multiplicity of charged particles in a collision event. Assuming that the number 
of particles hitting a silicon pad follows a Poisson distribution, the number of recorded hits 
(Nhn) in the r\ bin can be fitted as a function of the sum of the normalized energy deposits(-E) 
Nhu = N max {l — e~ E ' Emax ) to extract two parameters N max and E max . The parameter N max 
can be regarded as the effective number of active pads in the investigated rj bin and E max /N max 
is then taken as the average energy deposited by one charged particle in this rj bin. Thus the 
multiplicity in the investigated bin can be estimated as N = E / (E max / N max ) . From estimated 
charge multiplicities in two symmetric bins, N\ and N2, the variable C is calculated for each 
event. The rj bins are chosen to have equal size and to be symmetric about r/=0 in order to 
assure similar multiplicities in forward and backward regions. The RMS of the C distribution in 
each centrality and pseudo-rapidity bin is calculated directly as the fluctuations variable cr(C). 

Some preliminary <r(C) results^ are shown in Figure 1 and 2. Figure 1 shows the <r(C) 
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Figure 2. c(C) versus Ar; with r? bin centered at 2.0. All values of cr(C) were calculated using 
50% azimuthal acceptance. The estimated systematic error of 5% is not shown in the plot. 



dependence on the separation in ij space with fixed rj bin size 0.5. The data and reconstructed 
simulation (HIJING|2]+GEANTjS]) are similar. Both have a non-monotonic dependence on rj 
because of the presence of the detector effects. The HIJING cr(C) are slightly greater than 1, 
with the difference between the c(C) and 1 increasing with the ij separation. This indicates that 
the dynamic fluctuations of HIJING are greater than 0. The suppression in the mid-rapidity 
region is due to the fact that negative long range effect of the short range correlated particles 
selected into different bins cancels the positive short range correlations effect. Figure 2 shows 
the cr(C) dependence on the rj bin size with the r] centered at 2.0. The data and reconstructed 
simulation(HIJING+GEANT) are again similar. The HIJING <r(C) increases with the rj bin 
size. Such dependence on the bin size is a typical feature of particle correlations with finite 
correlation length. The data and reconstructed simulation c(C) also have the same trend of 
increasing with r/ bin size. 

Since detector effects dominate in the reconstructed multiplicity fluctuations as shown in 
Figure 1 and 2, it is necessary to clearly understand such effects in order to extract the dynamic 
fluctuations from the experimental measurement. Extensive GEANT detector simulations with 
different settings were used to investigate the relative importance of various sources in the 
observed detector effects. Four main detector effects were identified. The acceptance gap effect 
is significant in mid-rapidity region, but negligible in higher rj region. Such acceptance gap effect 
in <j{C) was caused by the asymmetric acceptance of the dead channels and gaps in eta coupled 
with the finite extension of the event vertex position. Since the energy deposit was used in 
the evaluation of the charged particle multiplicity, dE/dx fluctuations also contribute to a^ et - 
The dE/dx variations are due to Landau fluctuations and the wide range of detected particle 
velocities {[5 variations). The (3 variation contribution is more significant in the mid-rapidity 
region than in the high rj region. The Landau fluctuations contribution is nearly constant 



throughout r\ range covered by the octagonal detector. The contribution of secondary particles 
originating from primary particle decays and interactions is important primarily in the high 7] 
region because of the longer path length and higher probability of decay or interaction. 

To extract reliable dynamic fluctuations from the data, the analysis procedure was refined 
to suppress the detector effects as much as possible. For those residual detector effects which 
can not be suppressed through the improved analysis procedure, the simulations were used to 
reproduce them. In the mid-rapidity region and central collisions, the acceptance gap effect 
is the dominant detector effect. Through careful analysis of the cause of such detector effect, 
we designed an effective procedure for the suppression of the acceptance gap-induced detector 
effects by subtracting an event-by-event C offset, which is the mean of the C distribution in 
each (Z vtx , rj, Centrality bin). This C offset does not change the fluctuations due to other 
sources. Another improvement in the analysis is to use 77-dependent dE/dx hit cuts to suppress 
detector effects from the secondary contribution in high eta region. Using the improved analysis 
procedure, the detector effects in a(C) are greatly reduced and significantly less than those 
shown in preliminary results in Figure 1 and Figure 2. 

Three different simulations were used to evaluate the detector effects for the extraction of 
dynamic fluctuations: standard HIJING+GEANT, standard AMPT (7j +GEANT and modified 
HIJING+GEANT. The modified HIJING uses the HIJING event generator to generate events 
and randomizes the sign in 77 of the particles, thus destroying correlations between the particles 
while keeping other physics properties unchanged. Equation (2) was used to extract the 
detector effects from these simulations. The extracted detector effects are similar but have some 
small model dependence. Careful examination of the extracted detector effects' dependence 
on the model dynamic fluctuations helped identifying a linear relation between the extracted 
detector effects and the dynamic fluctuations of the model used in detector simulation, which is 
parametrized as 

°det = Vdetoi 1 - a(J dyn), ( 4 ) 

where the two parameters o"^ ef0 and a were estimated by fitting the a\ et as a function of o\ 
for each 77 bin. Such function fit was possible since there were three different simulations. The 
fitted (J 2 deia and a were then used to construct adjusted detector effects for the extraction of 
dynamic fluctuations with the following equation: 

cr 2 (C) = a 2 stat + a\ yn + (T 2 dem (l - aaj yn ). (5) 

Since the multiplicity was calculated using only half of the azimuthal acceptance, the 
extracted dynamic fluctuations were multiplied by the factor 2 to extrapolate dynamic 
fluctuations in half acceptance to that corresponding to full azimuthal acceptance. Such 
extrapolation was validated by the analysis of the dynamic fluctuations in simulations. In 
both HIJING and AMPT, the ratio of dynamic fluctuations for full azimuthal acceptance and 
dynamic fluctuations for half azimuthal acceptance in all the investigated 77 bins is 2. 

The dynamic fluctuations extraction procedure was applied to the reconstructed simulations 
to extract their dynamic fluctuations. As shown in Figure 3, the extracted dynamic fluctuations 
are consistent with the true dynamic fluctuations in all three simulations. Thus the analysis 
procedure for the extraction of dynamic fluctuations from reconstructed multiplicity was 
validated by simulations. We are now applying this procedure to extract the dynamic 
fluctuations from the Au+Au data. The analysis results compared with model predictions 
will be available in a later paper. We expect the new dynamic fluctuations results will provide 
useful insights on the intrinsic particle correlations in Au+Au collisions. 

In summary, we analyzed the multiplicity fluctuations in the inclusive charged particle 
production in Au+Au collisions at ^/snn = 200GeV within the pseudo-rapidity range of 
— 3 < 77 < 3. The preliminary fluctuations in the data are similar to that in reconstructed 
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Figure 3. Comparison of extracted and true model dynamic fluctuations in 0-20central 
collisions. The dynamic fluctuations analysis procedure as described in the text was used to 
extract the dynamic fluctuations from the reconstructed multiplicity of the simulations. 



simulations(HIJING+GEANT). The HIJING multiplicity fluctuations grow with increasing rj 
bin separation and width, revealing the existence of particle correlations in Au+Au collisions. 
An improved analysis procedure for extracting the dynamic fluctuations from the data was 
developed. This procedure was validated with simulations and it is now being used to extract 
the dynamic fluctuations from experimental data. The new data results will be available in 
near future. We expect these results will provide helpful insights on the intrinsic correlations in 
particle production. They could also be very useful input for theoretical model development for 
understating the particle production mechanism in heavy ion collisions. 
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